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ABSTRACT 

Chlorine and molecular hydrogen are known to be tightly linked together in the cold phase of the local interstellar medium through 
rapid chemical reactions. We present here the first systematic study of this relation at high redshifts using EE-bearing damped Lya 
systems (DLAs) detected along quasar lines of sight. Using high-resolution spectroscopic data from VLT/UVES and Keck/HIRES, 
we report the detection of Cli in 9 DLAs (including 5 new detections) out of 18 high-; DLAs with (V(FL) > 10 173 cnr 2 (including 
a new EL detection at z = 3.09145 towards J2100-0641) and present upper limits for the remaining 9 systems. We find a ~5cr 
correlation between A(Cli) and IV(EL) with only ~0.2 dex dispersion over the range 18.1 < log (V(EL) < 20.1, thus probing column 
densities 10 times lower those seen towards nearby stars, roughly following the relation (V(Cli) a 1.5 x ICE 6 x AffEL). This relation 
between column densities is surprisingly the same at low and high redshift suggesting that the physical and chemical conditions are 
similar for a given EL (or Cl i) column density. In turn, the AffCliJ/AffEL) ratio is found to be uncorrelated with the overall metallicity 
in the DLA. Our results confirm that neutral chlorine is an excellent tracer of molecule-rich gas and show that the molecular fraction 
or/and metallicity in the EL-bearing component of DLA could possibly be much higher than the line-of-sight average values usually 
measured in DLAs. 
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1. Introduction 

It has been shown that the global star-formation rate in the 
Universe gradually increases from z~10 to z~2-3 and then 
steeply decreases till the present epoch, z = 0 (see e.g. lDunlopl 
12011 . and references therein). Because metals are produced by 
stars, the determination of metal abundance in the gas pro- 
vides comple mentary information about star-formation history 
dRafelski et al.ll2012i) . This can be done using Damped Lyman- 
a absorption systems (DLAs) that represent the main reser¬ 
voir of neutral gas at h igh redshift dProchaska & Wolfel 12009: 
iNoterdaeme et al.l l2009) and are li kely to be located in g alaxies 
or in their close environment (e.g. lKrogager et ah 2012^. DLAs 


arise mostly in the warm neutral medium (e.g. iPetitiean et al.1 
2000; Ka nekar et al.ll2014) and have a multicomponent velocity 
structure, with metaljibsorption lines spread typically over 100- 
500 kms -1 (iLedoux et al.lIl998l> . In a small fraction of DLAs, 
the line of s ight intercepts cold gas, as traced by molecular hy- 
drogen (e.g. iNoterdaeme et al.l2008 .201_1: Balashev et al.120141) 
and/or 21-cm absorption (e.g. ISrianand et al.l 12012 ). Important 
progress has been made towards understanding the properties of 
the gas (through e.g. deriving physical conditions | Srianand et al.l 
2005; Noterdaeme et al.||2007at Oorgenson et al.ll2009 and phys¬ 
ical extent lBalashev et al.ll201 II) and the incidence of cold gas in 
DLA s has been related t o other properties (s uch as the metallic- 
ityj Petitiean et alj2006b or the dust content (iLedoux et al.l2003l 


iNoterdaeme et al.ll2008l) ). However, due to the strong saturation 
of H i Lyman series lines, it remains impossible to directly deter¬ 


mine the H i column density associated with the individual cold 
gas components traced by H 2 absorption. Even for metals, whose 
absorption lines are not saturated, it is very difficult to determine 
what fraction originates from the cold phase. Difficulties arise as 
well with the 21-cm absorption that do not always exactly coin¬ 
cide with H 2 absorption (iSrianand et al.|[20~i~3i ) although it could 
be due to the different structures of the optical and radio emitting 
regions of the background quasars. Out of all the metals, chlo¬ 
rine shows a unique behavior in the presence of H 2 . Because the 
ionization potential of chlorine (12.97 eV) is less than that of 
atomic hydrogen, chlorine is easily ionized in the diffuse neu¬ 
tral medium. However, this species reacts exothermically with 
H 2 at a very high rate converting rapidly Cl + into HC1 + . The lat¬ 
ter s ubsequently releases neutral chlorine through several chan¬ 
nels (iJuralll974clNeufeld & Wolfireil2009l) . This process is so ef¬ 
ficient that chlorine is completely neutral in presence of a small 
amount of H 2 . In our Galaxy the fact that chlorine abundance 
anti-c orrelates with the average number density al ong the line of 
sight dHarris & Bromagell 19841 Llenkins et all 19861 has been in¬ 
terpreted as chlorine depletion. However, models predict as well 
as observations indicate, that gas with moderate d ust c ontent 
presents negligible depletion of chlorine (e.g.|Neufeld & Wo I fire] 

l2009h ISavage & Sembachll996tlJenkinsl2009l) . Observationallv. 

a tight rel ation is indeed found between Cl 1 and Ho in the lo¬ 
cal IS M (l.lural 11974 ISonnentrucker et all 12004 iMoomev et al.l 
12012|) . In this letter, we present the first study of this relation at 
high redshift and over a wide range of column densities. 
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2. Data sample and measurements 


Since the first detection by iLevshakov & Varshalovichl (119851) . 
about two dozen H 2 absorption systems have been de¬ 
tected at high redshifts in QSO spectra. The detection 
limit of the stro ngest Cl 1 absorption line (1347A, f=0.0153 
dSchectman et alJFl993 111 in high quality spectrum (S/N ~ 50, 
R ~ 50000) corresponds to IV(Cli) ~ 10 12 c m~ 2 , The solar 
abundance of chlorine is 10" 6 ' 5 that of hydrogen (Asnlund et al.l 
2009) and given prev ious measurements of N(C1 i)/ 2N(H2) (e.g. 


Moomev et al.lf20T2h we conservatively limit our study to sys 
terns with MH 2 ) > 10 17 '"" 2 
Redshifts, H 1 and 


cm 

H 2 column densities, and metallici- 
ties were mainly taken from the literature and are based on 
VLT/UVES, Keck/HIRES or HST/STIS data. We refitted H 2 
absorption systems towards Q 2123-0050 and QJ2340-0053 
to take into account the positions of the detected Cl 1 com¬ 
ponents. We also detect a new H 2 absorption sys tem in the 
Z = 3.09 DLA towards J 2100-0641 in which Ijorgenson et al.l 
d2010l) have reported the presence of neutral carbon. C 1 is indeed 
known to be an excellent indicator of the presence of molecules 
(e.g. ISrianand et al.l 12005 ). We used the MAKEE package (T. 
Buries) to reduce archival data from this quasar obtained in 
2005, 2006 and 2007 under programs U17H (PI: Prochaska), 
G400H (PI: Ellison) and U149Hr (PI: Wolfe). We have found 
strong H 2 absorption lines from rotational levels up to / = 5 
(see Fig. Q} at z = 3.091485 with a total column density of 
log A(H 2 ) = 18.76 + 0.03. 

For all systems we retrieved data from the VLT/UVES 
or the Keck/HIRES archives. We reduced the data and fitted 
the lines using profile fitting. Neutral chlorine is detected in 
nine DLAs (Fig [2}. Four detections were already reported in 
the literature: Q 1232+ 082 dBalashev et al.ll2 01 ll). O 0812-3208 
dProchaska et al.ll2003l) . Q 1237+0647 ( Noterdaeme e t al. 2010) 
and Q 2140-0321 (Noterdaeme et al. submitted). The remaining 
five are new detections. We measured upper-limits of N(C1 1 ) for 
the remaining nine systems. We used mainly the 1347A Cl 1 line. 
Whenever possible, we also used Cli lines at 1088A, 1188A, 
1084 A, 1094A and 1085A, wi t h oscillator stre n gth from respec ¬ 
tively Schectman et al. (Il993l). lMortqn|j 2003|uMortqn| j2003l) . 
ISonn entruck er et al.l (l2006l) ~and lOliveira & Hebrardl ( 20061) . 

Table Q] summarizes the results of Cli measurements. We 
have kept all components with log MH 2 ) > 17. We did not use 
two known H 2 absorption systems towards Q 0013-0029 and 
J 091826.16+163609.0 since H 2 column densities in these sys¬ 
tems are not well defined. 


3. Results 

Fig. [3] shows the Cl 1 column density, MCI 1 ), versus MH 2 ) and 
compares our high-z measurements to those obtained in the local 
ISM using the Copernicus satellite dMoomev et al.l20l2l ). As can 
be seen our high-z measurements extend the relation to H 2 col¬ 
umn densities ten times smaller than those measured in the local 
ISM. Cl 1 and H 2 are found to be very well correlated (r = 0.95) 
over the entire MH 2 ) range. It is striking that measurements at 
high and low redshifts are indistinguishable in the overlapping 
regime (logMH 2 ) ~ 19-20.2). The correlation is seen over about 
three orders of magnitude in column density with a dispersion 
of about 0.2 dex only. A least-squares bisector linear fit provides 
a slope of 0.83 and 0.87 for the high-z and z = 0 data, respec¬ 
tively, with an almost equal normalization (logMCI 1 ) « 13.7 
at log MH 2 ) = 20). We note that the upper limits on MCI 1 ) lie 
mostly at the low MH 2 ) end and are least constraining since they 
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Fig. 1 . Voigt profile fits to the newly detected H 2 absorption lines from 
rotational levels J=0 to J=5 at z = 3.091485 towards J 2100-0641. The 
column densities are indicated (in logfcm 2 )) in the top right corner of 
each panel. 
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Fig. 2. Voigt profile fits to Cli 41347 absorption lines associated with 
high-z strong H 2 absorption systems. 


are compatible with the values expected from the above relation. 
For this reason, we will not consider them further in the discus¬ 
sion but still include them in the figures for completeness. The 
slopes are less than one, meaning that the Cl 1 /H 2 ratio slightly 
decreases with increasing MH 2 ). This is unlikely to be due to 
conversion o f Cl into H?C1 + and/or H C1, since chlorine chem¬ 
istry models (Neuf eld & Wolf ire 2009) as well as measurements 


(e.g. towards Sgr B2(S), Lis et al. l2010i) show that in diffuse 
molecular clouds only ~ 1% of chlorine is in the molecular form. 
A < 1 slope could in principle be due to dust depletion. However, 
there is no trend for increasing Cl depletion with increasing H 2 
or Cl 1 column densities in Galactic clouds (see IMoomev et al] 
l2012h . In addition, for high redshift measu rements, elemental 
abundance pattern (j Noterdaeme et al.l 2008|) as well as direct 
. iNotei 


measurements (e.g. iNoterdaeme et alJ |201 0|) indicate A„ < 0.2 
when modeling of chlorine chemistry (iNeufeld & Wolfire;l2009t) 
shows that for such low extinction (A v < 1) all chlorine is in the 
gas phase. 

A possibility is thus that the molecular fraction in the gas 
probed by Cl 1 is slightly increasing with increasing MH 2 ). It can 
be expected since H 2 self-shielding increase while Cl 1 is already 
completely in the neutral form. Finally, the similarity between 
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Table 1. Measurements of Cl i in strong Hi absorption systems at high redshift. 


Quasar 

£em 

£dla 

logWHi) [X/H] 

X 

Ref. 

zm 

log WHi) 

log WC1 1 ) 

b (kms *) 

[Cli/Hi] 

Q 0027-1836 

2.56 

2.40 

21.75+0.10 -1.63+0.10 

Zn 

1 

2.40183 

17.30+0.07 

< 12.71 

— 

< 1.61 

Q 0405-4418 

3.02 

2.59 

21.75+0.10 -1.12+0.10 

Zn 

2 

2.59475 

18.14+0.07 

< 12.71 

— 

< 0.77 

Q 0528-2505 

2.77 

2.81 

21.35+0.07 -0.91+0.07 

Zn 

3 

2.81098 

18.11+0.02 

11.92+0.08 

4.1 + 1.5 

0.01+0.08 







2.81112 

17.85+0.02 

11.73+0.11 

4.2+2.0 

0.08+0.11 

Q 0551-3637 

2.32 

1.96 

20.70+0.08 -0.35+0.08 

Zn 

4 

1.96214 

17 42 +0 - 45 

1 / Q 73 

< 12.40 

— 

< 1.54 

QJ0643-5041 

3.09 

2.66 

21.03+0.08 -0.91+0.09 

Zn 

5 

2.65860 

18.54+0.01 

12.51+0.05 

5.8+1.4 

0.17+0.05 

QJ0812+3208 

2.7 

2.63 

21.35+0.10 -0.81+0.10 

Zn 

6 

2.62628 

18.84+0.06 

12.79+0.05 

2.0+0.6 

0.15+0.08 







2.62644 

19.93+0.01 

13.78+0.27 

0.17+0.05 

0.05+0.27 

QJ0816+1446 

3.84 

3.29 

22.00+0.10 -1.10+0.10 

Zn 

7 

3.2874 

18.62+0.18 

< 13.65 

— 

< 1.23 







3.28667 

17.60+0.27 

< 12.76 

— 

< 1.36 

Q 1232+0815 

2.57 

2.34 

20.90+0.08 -1.35+0.12 

s 

8 

2.33772 

19.57+0.10 

13.49+0.08 

0.8+0.2 

0.12+0.13 

QJ1237+0647 

2.78 

2.69 

20.00+0.15 +0.34+0.12 

Zn 

9 

2.68959 

19.20+0.13 

13.01+0.02 

4.5+0.4 

0.01+0.13 

Q 1331+170 

2.08 

1.78 

21.18+0.04 -1.22+0.10 

Zn 

10,11 

1.77636 

19.71+0.10 

12.87+0.02 

5.7+0.5 

-0.64+0.10 

Q J1439+1118 

2.58 

2.42 

20.10+0.10 +0.16+0.11 

Zn 

12 

2.4184 

19.38+0.04 

< 13.25 

— 

< 0.07 

Q 1441+2737 

4.42 

4.22 

20.95+0.08 -0.63+0.10 

s 

13 

4.22401 

18.05+0.05 

< 12.86 

— 

< 1.01 







4.22371 

17.91+0.03 

< 12.66 

— 

< 0.95 

Q 1444+0126 

2.21 

2.09 

20.25+0.07 -0.80+0.09 

Zn 

14 

2.08696 

18.16+0.11 

< 12.42 

— 

< 0.46 

QJ2100-0641 

3.14 

3.09 

21.05+0.15 -0.73+0.15 

Si 

15 

3.09149 

18.76+0.03 

< 12.86 

— 

< 0.3 

QJ2123-0050 

2.26 

2.06 

19.18+0.15 -0.19+0.10 

S 

16 

2.05933 

18.09+0.02 

12.27+0.06 

2.6+0.5 

0.38+0.06 

QJ2140-0321 

2.48 

2.34 

22.40+0.10 -1.05+0.13 

P 

17 

2.33995 

20.13+0.07 

13.67+0.15 

5-10 

-0.26+0.18 

QJ2340-0053 

2.09 

2.05 

20.35+0.15 -0.92+0.03 

Zn 

15 

2.05456 

18.07+0.06 

12.25+0.15 

0.8+0.6 

0.38+0.16 







2.05473 

18.25+0.05 

12.13+0.18 

1 . 0 + 1.0 

0.08+0.19 

Q 2348-0108 

3.01 

2.43 

20.50+0.10 -0.62+0.10 

s 

18,19 

2.42688 

18.12+0.37 

< 13.86 

— 

< 1.94 







2.42449 

17.52+0.80 

< 13.13 

— 

< 0.81 


References. (1) INoterd aeme et al l d2007 al): (2) ILedoux et al.1 d 2003l): (3) INoter daeme et~aT1 J20 08I) : (4 ) ILedoux et al. I ]20(L g _ (5) 


lAlbOTiozVas^uezetal. 1 20141) : 161 Jorgensonet aljj|2009h: (7) Guhnaraeset al]j2012h: (81 Balashev^ al]j 201 lh: (9) NoterdaemeetaL 1 2010) 


(10) ICarswell et al.1 d201 1 iT T Tl l) iBaSev et alj <t20l71i ~Tl 2) ISrianand et al.1 J2008I) : ( 13) ILedoux et all <20061) : (1 4) ILedoux e t a ll (120031): ( 15) 
Jorge nson et all (l2010h : ( lfh lMalec et al.1 d20 10m (17) Noterdaeme et al. submitted; (18) lPetitiean et all d200fil) : (19) INoterdaeme et al. (l2007bl) . 



Fig. 3. Column densities of Cli versus that of H 2 . The red and blue 
points indicate the measurements at high redshift (this work) and in 
our Galaxy (tMoomev et al J l2C)T2h . respectively. The straight (dashed, 
dotted) lines show the respective least-squares bisector fits to the data. 


our Galaxy and high-z measurements at log WH 2 ) >19 might 
indicate that the chemical and physical conditions in the cold gas 
can be similar, otherwise fine tuning would be required between 


the different factors that impact on the N(C\ i) to WHi) ratio (e.g. 
number density, metallicity, dust content and UV flux). 

Before continuing further, we note that in Hz-bearing gas, 
chlorine is found exclusively in the neutral form (i.e. MCI) = 
(V(Cli)) (e.g. Q uralll97 4). Since we expect that all chlorine is in 
gas-phaseQ, the abundance of chlorine, [Cl/H], in Hi-bearing gas 
can be expressed as 

[Cl/H] = [C1 i/H 2 ] +log/, (1) 

where 

and / = 2WH 2 )/(2WH 2 ) + WHi)) is the molecular frac¬ 
tion. Therefore the ratio [C1 i/H 2 ] gives a direct constraint on 
the chlorine-based metallicity of H 2 -bearing gas provided the 
molecular fraction of H 2 -bearing gas is known. Conversely, if a 
constraint can be put on the actual chlorine abundance, [Cl i/H 2 ] 
can provide an estimate of the amount of Hi present in Hi- 
bearing gas. If Cl is depleted onto dust grains then the mentioned 
estimates of metallicity and molecular fraction will have to be 
corrected from the Cl depletion factor. 

The Fig.[4]shows [Cl i/H 2 ] as a function of the overall metal¬ 
licity for DLAs (given in Table Q} at high redshift or as a func¬ 
tion of [Cl/H] for clouds in our Galaxy. Since / < 1, [C1 i/H 2 ] 
gives an upper limit on the metallicity in H 2 -bearing gas, which 
is found to be roughly equal or less than solar metallicity. For 

1 We note that the presence of Cl n in the outer envelope of the H 2 
cloud is not excluded, but it does not influence our derivation 
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[X/H] dla , [C1 /H] mw 

Fig. 4. [Cl 1 /H 2 ] as a function of the overall metallicity for high-z DLAs 
(red points) and the chlorine-based metallicity for Milky-Way clouds 
(blue points). The dashed line represents the one-to-one relation. 

13 out of 21 Cl 1 bearing clouds in our Galaxy associated Cl n 
was measured (IMoomev et al.i 12012i) . Therefore we have esti¬ 
mated the overall chlorine abundance [Cl/H] of these clouds as 
(N(C 1 1 ) + tV(Cl ii))/(tV(Hi) + 2/V(H 2 )) (shown by blue circles in 
Fig. a. Unfortunately, for high redshift DLAs, not only Cl 11 is 
not detected but also DLAs contain several Hi clouds so that 
chlorine abundance of the very cloud of interest cannot be mea¬ 
sured. Therefore we consider the overall metallicity (averaged 
over velocity components) measured using another non-depleted 
element (usually Zn or S, see Table U}. In Fig. |4] it can be seen 
that the [Cl i/Ho] ratio is likely not correlated with the overall 
metallicity of the DLA (Pearson correlation coefficient 0.3 at 0.3 
significance level). For Milky-Way clouds it can be seen that 
[Cl/H] is typically one third solar, which can be interpreted as 
evidence for chlorine depletion ( Moomev et al]|20~i~2l) . 

The large difference between [Cl i/fL | and [X/H]dla for the 
high redshift clouds can be explained either by a molecular frac¬ 
tion / < 1 in H 2 -bearing clouds or by a higher metallicity in the 
Hi-bearing gas compared to the overall DLA metallicity or by 
both effects. If we assume that the metallicity in the Hi-bearing 
gas is equal to the overall DLA metallicity we find (using Eq.[T} 
that the molecular fraction in the Hi-bearing gas is typically an 
order of magnitude higher than the overall inferred DLA molec¬ 
ular fraction. Interestingly, two systems sitting close to the one- 
to-one relation are those where CO molecules have been de¬ 
tected (Q 1439+1118 and Q 1237+0647). In such systems, the 
H 2 component is probably fully molecularized and its metallic¬ 
ity is close to the overall metallicity of the DLA. 

4. Conclusion 

We have studied the neutral chlorine abundance in high red¬ 
shift (z~2-4) strong H 2 bearing DLAs with log ATH 2 ) > 17.3. 
These systems arise in the cold neutral medium of galaxies in the 
early Universe. We have used 17 systems from the literature and 
also present a new H 2 detection at z = 3.09145 in the spectrum 
of J 2100-0641. We have detected Cli absorption lines in half 
of these systems, including 5 new detections. The derived up¬ 
per limits for N(C1 1 ) for the remaining systems are shown to be 
consistent with the behavior of the overall population. Our mea¬ 
surements extend the Cl 1 -H 2 relation to lower column densities 
than measurements towards nearby stars. We show that there is 


a 5 cr correlation between the column densities of both species 
over the range 18.1 < log A'HU) < 20.1 with indistinguishable 
behavior between high and zero redshift systems. This suggests 
that at a given ATHi) the physical conditions are likely similar 
in our Galaxy and high-z gas, in spite of possible differences in 
the dust depletion levels. As we expect the Cl to be depleted less 
in the high-z absorbers studied here, we use the abundance of 
chlorine with respect to H 2 to constrain the molecular fraction 
and the metallicity in H 2 -bearing gas. Our results suggest that 
the molecular fraction and/or the metallicity in the H 2 and Cl 1 
bearing components could be much higher than the mean value 
measured over the whole DLA system. This implies that a large 
fraction of Hi is unrelated to the cold phase traced by H 2 . Fi¬ 
nally, our understanding of the formation of H 2 onto dust grains, 
self-shielding and lifetime of cold diffuse gas would certainly 
benefit from further observations of chlorine and molecular hy¬ 
drogen in different environments and over a wide range of col¬ 
umn densities. 
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